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Surface imaging of HD 199178 (V1794 Cygni) 
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Abstract. We present surface temperature maps for the FK Comae-type star HD 199178 (V1794 Cygni) calculated 
from high resolution spectra obtained in 1994 and 1995. The spot pattern evolves, but all maps reveal a large cool 
spot remaining nearly at the same high latitude. The main spot is 1200 - 1600 K cooler than the mean surface 
temperature. The observed slightly flat bottomed absorption lines would usually be interpreted as evidence for a 
large cool polar spot. We argue that antisolar surface differential rotation offers a better explanation for the box 
like shape of the line profiles. However, we do not find conclusive evidence for antisolar differential rotation and 
note that there are still other possible explanations for the slightly flat bottomed line profiles. 
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1. Introduction 

Doppler or surface imaging is, at present, the most power- 
ful tool to map spots on rapidly rotating late-type stars. 
Even so, the Doppler imaging maps should always be stud- 
ied critically. The method suffers from difficulties in mod- 
eling spectral lines in late-type stars, as well as from uncer- 
tainties in stellar and spectral parameters. These problems 
have been amply illustrated by conflicting results. For in- 
stance, the temperature maps obtained with some Doppler 
imaging methods most often show very large cool spots 



clear variations in their photometric rotation periods 



coverin g the p oles (cf. Vogt 1988; Strassmeier |1 994[ Hatzes 
& Vogt 1992| ), whil e the results with the m ethod s devel- 
oped by Piskunov ( 1991 ) and Berdyu gina ( |l998| ) rarely 



(Strassmeier et al. 1989; Cutispoto 1998). Variable peri- 
ods have also been discovered in the Can H&K emission 
(Donahue et al. 1996| ) . These periodicity changes are usu- 
ally explained as a signature of surface differential rotation 
combined with latitudinal activity migration. Apart from 
distorting the spot pattern, differential rotation would al- 
ter the shape of the line profile in the same way as a large 
cool or hot polar spot. This emphasizes the necessity of 
implementing surface differential rotation into the Doppler 
imaging techniques. 

The generally accepted empirical result is that differ- 
ential rotation is weaker in more rapidly rotating stars 



show polar spots (cf. Pisku nov e t al. 
et al. 1991 ; Korhonen et al. 200C| ). Since partly the same 
stars, and once even the same observations (Strassmeier 
et al. 1991), have been modeled, it is evident that the re- 
sults depend either on the applied methods or the selected 
stellar and spectral parameters. 

There are several studies on how line and stellar pa- 
rameters influence Doppler imaging (cf. Unruh & Collier 
Cameron 1995Q . But the influence of one important phe- 
nomenon, namely surface differential rotation, has usu- 
ally been neglected. Many active late-type stars show 



1990| - Strassmeier (Hall 1991). Hence the effects of differential rotation have 
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usually been bypassed in Doppler imaging. Although an- 
tisolar differential rotation has previously been noted to 
alter the line profiles in a sim ilar way as a large cool po- 
lar spot, Hatzes et al. ( 1996 ) concluded that this would 
require an unreasonably strong differential rotation. 

HD 199178 (V179 4 Cyg ni) is a rapidly rotating single 
G5 subgiant (Herbig 1958 ), which was class ified as a FK 
Coma e- type star by Bopp & Rucinski ( 1981 ). Jetsu et al. 
(|1990D derived a photometric rotation period of P p hot = 
3"? 337484 ±0<?000043. Furthermore, they reported a 9.1 
year period in the mean brightness level and indications 
of a shorter cycle in the light curve amplitude. In a recent 
ana lysis o f photometry between 1975 and 1995, Jetsu et 
al. ( 19991: ) showed that the photometric period varied and 
interpreted this as a signature of differential rotation. 

We present new temperature maps of HD 199178 for 
the years 1994 and 1995. The maps were calculated with 
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the surface imaging technique developed by Piskunov 
(1991). Previous Doppler images of this star display a 
large cool polar spot and some smaller high and low lati- 
tude features (Vogt 1988; Strassmeier et al. 1999). Unlike 



i; Strassmeier et al. 
these earlier studies of HD 199178, we take into account 
the possible effect of differential rotation on the rota- 
tional broadening of photospheric absorption lines. We 
also demonstrate how different values for the differen- 
tial rotation coefficient alter the retrieved surface imaging 
map. 

2. Stellar parameters of HD 199178 

HD 199178 has an apparent visual magnitude of V « 7^2. 
The Hipparcos-satellite parallax measurement of tt « 
10.68 ± 0.73 mas gives a distance of 94 pc (Perryman 
1997 ). The long-term mean colour indices B— V ss 0.79 and 
V — R « 0.68 determined by Jetsu et al. ( 1990| ) suggested 
an average "background" stellar surface temperature of 
T G ff w 5375 K. The above results agre e with the G5iii-iv 
spectral classification by Herbig ( 1958 ). 

The new Hipparcos distance measurement allows us to 
estimate the absolute magnitude and radius of the star. 
We use 7^09 as the apparent visual magnitude of the 
unspotted star (Strassmeier et al. 1999 ). HD 199178 lies in 
the direction of the North America and Pelican Nebulae 
complex, but in the foreground. Assuming an interstellar 
absorption of 0'. n 2 (Straizys et al. |1989| : Fig. 2) yields an 
absolute visual magnitude of 2™02. A bolometric correc- 
tion of BC = -0.184 (Strassmeier et al. |1999| ) gives the 
bolometric magnitude Mboi ~ 1™84. Finally, if the unspot- 
ted temperature of T e ff = 5400 K is combined to the solar 
T = 5800 K and M bo i = 4 m 64, the Stefan-Boltzmann 
relation gives R « 4.2i? . 

The Barnes-Evans relation is based on the same prin- 
ciple as the above, except that the value of the bolometric 
correc tion is not required. We use the formulation by Lacy 
C|l977b : 



\g(R/R Q ) = 7.4724 - 0.2V - 2F V + Igd, 



(D 



where F v = 3.977 - 0.429(V - R) and [d] = pc. The 
combination V = ?. n 0, V - R = 0.68 and E V _ R = 0.05 
gives F v w 3.707 and R w (4.3 ± O.7)i? . (Here V is 
the mean apparent visual magnitude without interstellar 
absorption.) 

The rotation period P ps 3^3 and projected rotational 
velocity vsmi ps 70 km s _1 give a lower limit Rsini « 
4.6i? for the radius. There is thus a slight discrepancy 
between the different radius determinations. Anyway, the 
above estimates and the lower limit i?sinz suggest that 
4.6i?Q R ^ 5.Oi?0. This is equivalent to an inclination 
of 67° < i < 90°. The exceptional nature of HD 199178 
(rapid rotation, spots, possibly a coalesced binary) may 
mean that the equatorial radius will be underestimated. 
In our calculations we will thus allow for an inclination 
i > 50°. 

A variable photometric rotation period was dete cted in 
the long-term photometry by Jetsu et al. ( |l999b[ ). They 



analyzed 114 sets of photometric data obtained between 
1975 and 1995. Each set contained at least 8 nights of 
observations and had a length of about 30 nights (Jetsu 
et al. 1999a). The regular yearly change rate of the pho- 
tometric period was 3.3 %. In addition, irregular seasonal 
changes of 7.5 % were found. The varying photometric pe- 
riod was explained by surface differential rotation quanti- 
fied by |a| <; 0.075, where 



Q; (^equator ^polc)/^cquator- 



(2) 



Even this lower limit for the differential rotation is consid- 
erably larger than that expected from the results by Hall 
( 1991 ). If the latitudinal range of the spot activity does 
not extend all the way from the equator to the pole, then 
|a| could even reach the solar value of 0.2. We thus assume 
that the realistic differential rotation coefficient values are 
within the limits -0.20 < a < 0.20. 

Differential rotation would alter the shape of the line 
profile. The case of a > induces a "sharper" profile, 
while a < causes a more "flat bottomed" profile (cf. 
Bruning 1981). Thus a negative a modifies the line profile 
much in the same way as a cool polar spot. 

Since there is no unique photometric period, in prin- 
ciple different periods could be used for different sea- 
sons. But this would make it harder to compare the sur- 
face images. Therefore we chose to use the mean period 
P = 3"?3250 calculated for the years 1994 an 1995. Our 
ephemeris for the rotational phases was 



HJD min = (2449511.150 ±0.070) + (3.3250±0.0012).E, (3) 

where the HJD m j n is the epoch of the light curve minimum 
(Jetsu et al. |1999b[). 



3. The surface imaging technique 

Spots on the surface of a rapidly rotating star will cause 
distortions ("bumps") in the observed spectral line pro- 
files. As the star rotates, these bumps will move across the 
absorption line profiles. The main idea of surface imaging 
is to trace these distortions and retrieve a surface map 
of the star. We use the surface im aging technique devel- 
oped by Piskun ov (Piskunov et al. 1990 , Eqs. 3-6 and 8; 
Piskunov 1991, Eqs. 5-8), but also include Johnson V- 
photometry in the inversion. 

Local line profiles and F-band fluxes are calculated for 
a set of temperatures and limb angles using stellar model 
atmospheres. This table is then used for the disk integra- 
tion of a given surface temperature distribution. The sur- 
face imaging problem is solved by searching for a surface 
temperature distribution that minimizes the discprepancy 
between the observations and the calculated line profiles 
and light curve. 

Since we always have a confined temperature range 
of model atmospheres, the solution should be limited to 
temperatures between the minimum and maximum of the 
models (T m ; n and T max ). We do this by adding a "penalty 
function" to the quantity to be minimized. We use an 
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integral of a double-sided Heaviside step function, which is 
nearly constant between T m i n and T max , but grows nearly 
linearly when the temperature is below T m i n or exceeds 
Tmax- Thus the penalty function merely prevents extreme 
temperatures in the solution. 

The crucial point of surface imaging is how the line 
profiles for a given temperature distribution are calcu- 
lated. Errors in the calculated line profiles will intro- 
duce artifacts in the resulting map and prevent the solu- 
tion from reaching the observational signal-to-noise level. 
Errors in the line profiles can arise from numerous prob- 
lems, e.g. wrong line parameters, missing lines, non-LTE 
effects, wrong limb darkening or inadequacies of the model 
atmospheres. The fact that the lines used for Doppler 
imaging are strong absorption lines, whose cores are 
formed high in the photosphere, puts extra demands on 
the line modeling. 

The errors in the line profiles are not the only source 
of systematic errors. For example, an erroneous vsiiii will 
manifest its elf as a cool or hot belt in the image (cf. 
Berdyugina 1998). Errors in the microturbulence veloc- 



ity will mainly change the average temperature of the so- 
lution. Doppler images are not sensitive to errors in the 
macroturbulence velocity (cf. Unruh & Collier Cameron 
19951 ). A wrong inclination will cause a latitudinal dis- 
placement and a slight distortion of the spot. Combined 
with surface differential rotation, the inclination will also 
influence the line profile. Decreasing the inclination will 
increase the effect of differential rotation on the line pro- 



file (cf. Bruning |l98l| ). 

Differential rotation will, of course, also change the 
spot distribution, since different latitudes have different 
rotational periods. In this study we do not take this effect 
into account. Firstly, differential rotation becomes impor- 
tant in this sense only if the locations of spots extend over 
a large latitude range. Secondly, we assume that the mag- 
netic field strength is sufficient to resist the shearing effect 
of differential rotation within the spot(s). 

Since it is not (yet) possible to calculate perfect line 
profiles, one can always expect some artifacts in Doppler 
images. Fortunately, most of these artifacts can easily be 
identified, since certain shapes of artifacts are related to 
particular error sources. The same error in the overall 
shape of all calculated line profiles may manifest itself as a 
spot band at a certain latitude (or a polar spot). However, 
if the rotational phase coverage is inadequate, this artifact 
might not be axisymmetric. Minor errors in the calcula- 
tions (e.g. missing lines in blends) or the data reduction 
(e.g. remnants of cosmic spikes or uncertainty in the con- 
tinuum normalization) may cause thin high contrast fea- 
tures, like vertical stripes, arches or even ovals. For exam- 
ple, the remains of a cosmic spike in one spectrum may 
cause a thin high contrast artifact at a stellar longitude 
corresponding to the wavelength of the spike and rotation 
phase of the observation. This in turn might cause addi- 
tional hot and cold features to cancel out the effect of the 
artifact in the remaining observed phases. A missing line, 
on the other hand, may lead to an artifact repeated at 



each observed phase. Noise in the spectra might be inter- 
preted as a series of small "bumps" . Because such bumps 
will appear at arbitrary wavelengths, the resulting arti- 
facts cannot be latitudinally confined, and will therefore 
be seen as vertical stripes in the image. 

It is clear, that the retrieved longitude for the spots will 
usually be more reliable than the latitude. In some cases 
the longitude of the spot can be estimated even with a 
single spectrum, while the determination of the latitude 
requires a good phase coverage. Moreover, the latitude is 
sensitive to errors in the inclination and the general shape 
of the rotationally broadened line profile. The influence of 
all error sources is greatly reduced by using many spectral 
lines and simultaneous photometry in the inversion. 

Table 1. The spectral observations of HD 199178: t = 
HJD— 2440000 for the middle of the exposures, rotational 
phases calculated from the ephemeris of Eq. (||), S/N ra- 
tios and the resolution element W\ 

July 1994, Rozhen, 14^=0.20 A 



t 


Phase 


S/N 


t 


Phase 


S/N 


9552.341 


0.388 


244 


9559.479 


0.535 


212 


9552.457 


0.423 


268 


9560.303 


0.783 


171 


9552.576 


0.459 


214 


9560.474 


0.834 


320 


9553.309 


0.679 


219 


9561.350 


0.098 


162 


9554.559 


0.055 


251 


9562.353 


0.399 


286 


9555.421 


0.315 


251 


9562.537 


0.455 


206 


9556.506 


0.641 


223 


9563.352 


0.700 


236 


9557.537 


0.951 


248 


9564.393 


0.013 


316 


9558.346 


0.194 


276 


9565.381 


0.310 


263 


9558.561 


0.259 


178 








August 1994, NOT-SOFIN/ 2nd camera,^ =0.10 A 


t 


Phase 


S/N 


t 


Phase 


S/N 


9578.544 


0.269 


209 


9584.412 


0.034 


284 


9579.608 


0.589 


281 


9585.409 


0.334 


238 


9580.521 


0.863 


277 


9586.410 


0.635 


308 


9581.503 


0.159 


262 


9587.441 


0.945 


263 


9582.418 


0.434 


246 


9588.424 


0.240 


263 


9583.412 


0.733 


282 








November 1994, NOT-SOFIN/ 1st camera, W A = 


0.05 A 


t 


Phase 


S/N 


t 


Phase 


S/N 


9669.355 


0.580 


150 


9675.323 


0.375 


210 


9670.305 


0.866 


190 


9676.322 


0.676 


210 


9671.412 


0.199 


180 


9677.328 


0.978 


230 


9672.310 


0.469 


220 


9678.325 


0.278 


250 


9673.309 


0.770 


210 








July 1995, NOT-SOFIN/ 2nd camera 


,Wa=0.10 


A 


t 


Phase 


S/N 


t 


Phase 


S/N 


9910.599 


0.135 


295 


9917.580 


0.235 


301 


9912.573 


0.729 


454 


9918.588 


0.538 


321 


9913.598 


0.037 


332 


9919.634 


0.852 


295 


9914.576 


0.331 


376 


9920.623 


0.150 


319 


9915.616 


0.644 


325 


9921.532 


0.423 


362 


9916.628 


0.948 


281 
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Table 2. Atomic lines, lower excitation potential and adopted gf values 



Line 


A (A) 


Xiow (eV) 


ig(s/) 


Line 


A (A) 


Xiow (eV) 


lg(ff/) 


Line 


A (A) 


Xiow (cV) 


ig(s/) 


Si I 


6407.070 


5.871 


-2.200 


Smi 


6425.905 


0.282 


-1.815 


Til 


6435.315 


2.250 


-2.875 


Zn 


6407.216 


0.154 


-2.700 


Si i 


6426.291 


5.954 


-2.470 


Vi 


6435.490 


2.684 


-1.582 


Fen 


6407.251 


3.889 


-3.699 


Smil 


6426.628 


1.746 


-0.559 


Cei 


6436.399 


0.467 


-0.822 


Si i 


6407.291 


5.871 


-1.500 


Zn 


6426.706 


2.746 


-0.620 


Fei 


6436.407 


4.186 


-2.360 


Co i 


6407.491 


4.395 


-1.471 


Si i 


6426.911 


5.954 


-2.900 


Yi 


6437.169 


2.294 


-0.520 


Fei 


6407.643 


4.076 


-3.620 


Smil 


6428.355 


1.370 


-1.150 


Til 


6437.610 


3.407 


-1.889 


Fei 


6408.018 


3.686 


-1.018 


Ndn 


6428.645 


0.205 


-1.831 


Eull 


6437.640 


1.320 


0.427 


Si 


6408.157 


7.868 


-1.510 


Cai 


6428.815 


4.441 


-1.184 


Si i 


6437.703 


5.863 


-2.210 


Fei 


6408.332 


4.386 


-3.563 


Fei 


6429.071 


4.294 


-3.407 


Nil 


6437.992 


5.389 


-2.984 


Sri 


6408.459 


2.271 


0.500 


Pru 


6429.629 


1.615 


-0.386 


Vi 


6438.088 


2.684 


-2.070 


Gdi 


6408.545 


0.213 


-2.052 


Ndi 


6429.840 


1.337 


-0.530 


Cn 


6438.586 


3.890 


-2.571 


Si i 


6408.671 


5.984 


-3.090 


Nil 


6429.859 


4.167 


-2.104 


Fei 


6438.755 


4.435 


-2.329 


Vi 


6408.814 


2.616 


-1.992 


Co i 


6429.906 


2.137 


-2.460 


Cai 


6439.075 


2.526 


0.450 


Til 


6409.390 


3.570 


-1.631 


Cai 


6430.127 


3.910 


-1.907 


Til 


6439.326 


3.590 


-1.318 


Vi 


6410.416 


2.581 


-1.889 


Co i 


6430.290 


4.049 


-1.028 


Fei 


6439.554 


4.473 


-4.304 


Lai 


6410.975 


0.373 


-0.950 


Nbi 


6430.443 


0.740 


-1.230 


Til 


6439.705 


0.813 


-6.458 


Sci 


6410.984 


4.967 


-2.870 


Vi 


6430.472 


1.955 


-0.600 


Cei 


6439.964 


0.294 


-0.752 


Fei 


6411.107 


4.733 


-1.905 


Cai 


6430.793 


3.910 


-2.129 


Si i 


6440.566 


5.616 


-2.480 


Vi 


6411.276 


1.950 


-2.059 


Fei 


6430.846 


2.176 


-2.106 


Ni 


6440.938 


11.764 


-1.140 


Cri 


6411.580 


3.892 


-2.478 


Smil 


6431.006 


1.359 


-1.158 


Til 


6440.947 


3.590 


-1.810 


Fei 


6411.649 


3.654 


-0.645 


Cai 


6431.099 


3.910 


-2.606 


Mm 


6440.971 


3.772 


-0.938 


Co i 


6411.884 


2.542 


-2.228 


Vi 


6431.623 


1.950 


-0.287 


Fei 


6441.780 


5.033 


-1.928 


Til 


6412.192 


3.581 


-1.480 


Pru 


6431.805 


1.416 


-0.339 


Si I 


6442.777 


6.125 


-1.240 


Fei 


6412.200 


2.453 


-4.441 


Smil 


6431.978 


1.413 


-1.460 


Fen 


6442.955 


5.549 


-2.885 


Til 


6413.108 


0.048 


-5.006 


Nil 


6431.994 


3.542 


-2.500 


Fei 


7507.261 


4.415 


-3.482 


Sci 


6413.324 


0.021 


-2.677 


En 


6432.531 


1.309 


-0.919 


Fei 


7508.533 


4.988 


-2.318 


Gai 


6413.468 


3.073 


-0.300 


Ndi 


6432.650 


1.287 


-0.160 


Fei 


7508.623 


4.143 


-2.741 


Ci 


6413.550 


8.771 


-1.380 


Fen 


6432.680 


2.891 


-3.540 


Til 


7509.728 


1.749 


-3.030 


Pr II 


6413.678 


1.132 


-0.799 


Vi 


6433.162 


1.945 


-1.308 


Si I 


7510.785 


5.984 


-1.770 


Mm 


6413.945 


3.763 


-1.910 


Nbi 


6433.204 


0.657 


-1.480 


Fei 


7511.020 


4.178 


0.059 


Si i 


6414.188 


6.083 


-2.450 


Si i 


6433.337 


5.614 


-3.210 


Ndn 


7511.125 


1.773 


-0.871 


Nil 


6414.581 


4.154 


-1.220 


Si i 


6433.457 


5.964 


-1.590 


Fei 


7512.115 


2.279 


-5.195 


Nil 


6414.796 


4.165 


-2.463 


Fen 


6433.814 


6.219 


-2.371 


Fei 


7512.139 


4.143 


-2.947 


Si i 


6414.980 


5.871 


-1.100 


Zn 


6434.388 


2.798 


-0.090 


Cm 


7512.745 


4.756 


-2.630 


Fei 


6415.505 


4.435 


-3.912 


Cel 


6434.388 


0.028 


-1.437 


Ndn 


7513.736 


0.933 


-1.241 


Si 


6415.522 


7.870 


-1.360 


Yi 


6435.004 


0.066 


-0.820 


Fen 


7515.091 


5.823 


-2.362 


Ndn 


6425.779 


1.649 


-0.935 


Vi 


6435.158 


1.942 


-1.458 


Fen 


7515.831 


3.903 


-3.432 



4. Observations 

High resolution spectra of HD 199178 were measured in 
July 1994, August 1994, November 1994 and July 1995. 
The July 1994 data set was obtained with the 2 m Ritchey- 
Chretien telescope coude-spectrograph of the National 
Astronomical Observatory (Rozhen), Bulgaria. The rest 
of the spectroscopic observations were made using the 
SOFIN high resolution echelle spectrograph at the 2.5 m 
Nordic Optical Telescope (NOT), La Palma, Spain. The 
observed spectral regions were 6409-6414 A (August 1994 
and July 1995), 6427-6441 A (all seasons) and 7509-7514 
A (August 1994, November 1994 and July 1995). The typ- 
ical signal-to-noise ratio was between 200 and 300. The 
spectral observations are summarized in Table |j 

The spectral obse rvatio ns were reduced with the 3A 
software system (Ilyin [2000 ). The reduction included bias, 
cosmic ray, flat field and scattered light correction, wave- 
length calibration and normalization. Special care was 
taken during the last two steps, where errors can cause 
artifacts in the surface images. 



The calibration of the wavelength scale was obtained 
from Th-Ar comparison spectra. For the echelle spectra, 
we used all available spectral orders of the comparison 
frames to calculate a two-dimensional (kX) wavelength so- 
lution, taking into account distortions of the line positions 
(Ilyin |2000p . Small shifts between the stellar and compar- 
ison spectra (e.g. due to flexure of the spectrograph) were 
corrected using atmospheric lines. The spectra were re- 
binned to a common wavelength grid. 

The photospheric absorption lines of HD 199178 are 
broad and it is hard to determine the real continuum in 
the spectrum. Therefore the continuum normalization had 
to be done in two steps. A "quasi-continuum" consisting 
of points near to the real continuum level was used as a 
first approximation. The final normalization was done by 
utilizing a rotationally broadened synthetic spectrum. 

Simultaneous or nearly simultaneous Johnson V- 
photometry for July 1994, August 1994 and July 1995 was 
obtained with the 60 cm telescope at the Mount Maidanak 
Observatory and the APT Phoenix- 10 tele scope (subsets 
SET=108, 109 and 111, Jetsu et al. |l999a|). 
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5. Temperature maps 

The following three spectral regions were used for 
surface imaging: 6409.0 - 6413.5 A (strongest line 
Fei A641 1.649 A), 6427.5 - 6441.0 A (Fei A6430.846 
A; Fen A6432.680 A; Cai A6439.075 A) and 7509.0 - 
7513.5 A (Fei A7511.020 A). The local line profiles were 
calculated for 20 limb angles and temperatures 3500 - 
6500 K (with a st ep of 250 K) using a code written by 
Bcrdyugina (1991) and atmospheric models by Kurucz 
(1993) with solar element abundances. Spectral line pa- 
rameters were obtained from the Vienna Atomic Line 



Database (Kupka et al. 1999). All significant atomic lines 
and the most important molecular lines (CN and TiO) 
were included in the calculations. The final adjustment of 
the gf values was done after analyzing the wavelength de- 
pendence of the deviation between a synthetic spectrum of 
HD 199178 and the mean of the observations. Of the 111 
lines used in the calculations, the gf value of 36 lines was 
changed. The changes in the strongest lines, Fe I A6411.649 
A, Fei A6430.846 A, Fen A6432.680 A, Cai A6439.075 
A and Fei A7511.020 A, were -0.1 < Alg(gf) < 0.02. It 
should be emphasized that the adopted gf values may be 
biased by uncertainties in element abundances. A list of 
all atomic lines that we used, can be found in Table ||. 

Table 3. Stellar parameters used for the surface images 
in Figs. |l| - | 



Parameter 


Adopted value 


Gravity lg(.g) 


3.5 


Inclination i 


60° 


Rotation velocity v sin i 


70 kms" 1 


Differential rotation a 


-0.17 


Rotation period P ro t 


3.3250 


Micro turbulence n m icro 


1.4 km s" 1 


Macro turbulence ti macro 


4.0 km s" 1 


Element abundances 


solar 




The temperature maps were derived using the 
INVERS7-9PD code, originally written by N.E. Piskunov, 
with some changes made by T. Hackman. We first tested 
different values of rotational velocity (68 km s _1 < vsmi< 
74 km s _1 ), microturbulence (1.2 km s _1 < tWcro < 1.8 
km s^ 1 ), inclination (50° < i < 85°) and surface differen- 
tial rotation (—0.20 < a < 0.20). Differential rotation was 
implemented by adjusting the angular rotational velocity 
at each latitude b using the Maunder rotation formula 

0(6) = Oequator (1 - a sin 2 b) . (4) 

The best solution (the smallest deviation between the 
spectroscopic observations and the calculated line profiles) 
was obtained with the stellar parameters given in Table [|. 

Before combining all available spectral regions, sepa- 
rate images were calculated for each line blend in order 
to check for inconsistencies. The individual maps for each 
spectral region gave the same main features, but contained 



Fig. 1. The July 1994 image of HD 199178 with calculated 
and observed spectra and photometric V-curve. The mean 
deviation of the spectroscopic observations from the model 
was d = 0.56%. 



more artifacts (small high contrast features, especially at 
low latitudes) than the maps obtained from the combina- 
tion of all spectral regions. 

All surface images show a large cool spot remaining 
at approximately the same latitude (Figs. |l|-| and Table 
|]). The mean temperature of the visible surface is 5300 
K. The polar region is about 4900 K. The lower latitude 
features in the images are probably artifacts. In all maps 
there are symmetry effects typical for surface imaging of 
stars with relatively high inclination (i > 60°), i.e. "shad- 
ows" of the high latitude spot below the equator. The 
bright features, reaching a maximum of T c g m 5800 K, 
are probably also artifacts. 
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Fig. 2. The August 1994 image of HD 199178 with calcu- 
lated and observed spectra and photometric 1^-curve. The 
mean deviation of the spectroscopic observations from the 
model was d — 0.43%. 



The July 1994 image shows a large high latitude cool 
spot (3700 K) with a longitudinal appendage (3900 K). 
In the August 1994 image, the coolest part of the spot 
(3700 K) is more concentrated and the appendage (4200 
K) is weaker. In November 1994, the spot appears to be 
partly disintegrated. The coolest part of the spot (4100 
K) seems to be at the position of the previous appendage. 
No photometric data was available for the November 1994 
image, but there were clear changes in the light curves for 
Jul y - Se ptember 1994 (subsets SET=108-110, Jetsu et 



al. |l999b|) , which support the spot evolution. The July 
1995 image resembles the July 1994 image, but with an 
additional new appendage stretched "southwards" . The 
coolest part of the spot (3700 K) has regained its strength. 

The best surface imaging solution was achieved with 
the differential rotation coefficient a = —0.17. This im- 
plies strong antisolar differential rotation in HD 199178, 
i.e. that the pole would rotate faster than the equator. It 
seems that current models for global circulation in outer 
stellar convection zones predict only solar- like differential 
rotation (cf. Kitchatinov & Riidiger 1999 ). Therefore we 




phase=O.O0j^pMean = 5296K phase=0.2 & c -p M _ 



6000K 
5850K 
5700K 
5550K 
5400K 
5250K 
5100K 




Fig. 3. The November 1994 image of HD 199178 with 
calculated and observed spectra (no ^-photometry avail- 
able). The mean deviation of the spectroscopic observa- 
tions from the model was d = 0.50%. 

Table 4. Temperatures T spot , longitudes I and latitudes 
b of the coolest concentrated parts of the spot in each im- 
age, and photometric periods P for each set from Jetsu 
et al. 1999b| . (The photometric periods are listed only for 
comparison, each image was derived using the same pho- 
tometric period P = 3.3250, i.e. the mean period for the 
years 1994-1995.) At phase the longitude I = -360°<?(> is 
at the center of the visible stellar disk 



Image 


T spot (K) 


en 


b(°) 


P(d) 


July 1994 


3700 


-5 


67 


3.352 


August 1994 


3700 


-15 


61 


3.384 


November 1994 


4100 


-70 


65 




July 1995 


3700 


-39 


68 


3.2829 



also wanted to test the possibility of no differential ro- 
tation and solar-like differential rotation. Our Fig. |s| dis- 
plays the July 1995 surface imaging solutions obtained 
with a — (rigid body rotation) and a — 0.15 (solar- 
like differential rotation). The best fits with a = and 
a = 0.15 were achieved by increasing usini to 72 km 
s^ 1 and 74 km s , respectively. The mean deviation (d) 
of the calculated line profiles from the observations was 
slightly higher for a = and a — 0.15 than for a — —0.17, 
i.e. the antisolar case provided a better fit to the data 
(Table |). 

The images in Fig. || demonstrate the difficulties in sep- 
arating the effect of a polar spot from that of differential 
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Fig. 4. The July 1995 image of HD 199178 with calculated 
and observed spectra and photometric V^-curve. The mean 
deviation of the spectroscopic observations from the model 
was d = 0.37%. 



rotation. With no differential rotation the polar region is 
~ 800 K cooler than the mean surface temperature, while 
a = 0.15 gives a cool polar spot with T c g ~ 3900 K, much 
like the results by Strassmeier et al. ( 1999[ ). The same 
general effect was verified for all other seasonal images. 
Except for the polar cap, the images in Fig. || are very 
similar. The coolest part of the spot remains at approx- 
imately the same position. Similar artifacts can also be 
seen in all three images, in particular some bright features 
and smaller spots near the equator. 



6. Differential rotation or polar spot? 

The surface imaging solution with antisolar surface differ- 
ential rotation might raise the question of whether this is 
a real phenomenon or just a trick to eliminate a cool polar 
spot. Thus, what we need is a way to determine which one, 
the polar spot or differential rotation, gives a better solu- 
tion. Whe n anal yzing the T Tauri star Sz 68, Johns-Krull 
& Hatzes ( 1997 ) tackled this problem by calculating sep- 
arate temperature maps for different spectral lines. They 





iiOOOK 
5850K 
5700K 
5550K 
5400K 
5250K 
5100K 
4950K 



Vsin(i)= 74.0 km/s , i= 60, a = 0.15 




Fig. 5. The July 1995 image with differential rotation co- 
efficients a — —0.17 (same as Fig. |i|), 0.0 and 0.15. 

Table 5. The quality of the July 1995 surface imaging 
solution using different values of a: mean deviation d of 
the calculated line profiles from the spectroscopic obser- 
vations and linear correlation coefficient r between next- 
to-neighbouring residuals 





vsini 


a 


d 


r 


70 


km s _1 


-0.17 


0.373% 


0.507 


72 


km s _1 


0.00 


0.388% 


0.540 


74 


km s -1 


0.15 


0.410% 


0.572 



argued for solid body rotation and a polar spot, because 
this alternative gave a larger correlation between the in- 
dividual maps. 

We tested whether this approach works in our case. 
First we simulated synthetic spectroscopic "observations" 
using as input a stellar surface with a single high lat- 
itude spot, but no polar spot, and stellar parameters 
wsini=70 km s~\ a = -0.2, i = 60° (CASE i). To 
imitate real observations, we added normally distributed 
noise (a = 0.30%) to the synthetic observations. These 
synthetic spectra were then used to calculate temperature 
maps with two sets of stellar parameters: vsini =70 km 
s _1 , a = —0.2, i = 60° and i;siiu=72 km s _1 , a = 0.0, 
i = 60°. Maps were calculated both using all spectral re- 
gions and separately for each spectral region (6411 A, 6431 
A, 6439 A and 7511 A). With the correct stellar param- 
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eters, we retrieved a surface imaging solution consistent 
with the input stellar surface. As expected, neglecting the 
differential rotation, i.e. using a = 0.0, resulted in a solu- 
tion with a cool polar spot in addition to the high latitude 
spot. 

For both sets of stellar parameters we calculated the 
correlations between the maps retrieved for the different 
spectral lines. The test results showed that the correla- 
tion was higher when the differential rotation was wrongly 
neglected (Tables || and [?]). The reason for this is clear. 
Neglecting a substantial differential rotation, one always 
obtains the same artifact, a polar spot, regardless of the 
spectral line used for the image. The presence of the same 
polar spot in all images will naturally increase the corre- 
lation between the images. The correlation only measures 
the consistency between the calculated images. It will not 
reveal the inconsistency between the original data (no po- 
lar spot) and the solution (polar spot). 

Table 6. CASE i: Correlation between solutions for dif- 
ferent spectral lines. The solutions were retrieved neglect- 
ing differential rotation (a — 0), although the original 
data was simulated with a differential rotation coefficient 
a = -0.2 





6411 A 


6431 A 


6439 A 


7511 A 


6411 A 


r=1.000 








6431 A 


r=0.975 


r=1.000 






6439 A 


r=0.959 


r=0.967 


r=1.000 




7511 A 


r=0.961 


r=0.953 


r=0.970 


r=1.000 



Table 7. CASE i: Correlation between solutions for dif- 
ferent spectral lines. Both the solutions and the original 
data were simulated with the differential rotation coeffi- 
cient a = —0.2 





6411 A 


6431 A 


6439 A 


7511 A 


6411 A 


r=1.000 








6431 A 


r=0.899 


r=1.000 






6439 A 


r=0.918 


r=0.907 


r=1.000 




7511 A 


r=0.887 


r=0.913 


r=0.892 


r=1.000 



An alternative approach to check the reliability of 
Doppler imaging maps, is to estimate the systematic er- 
rors in the calculated line profiles. As opposed to noise, the 
systematic errors of neighbouring wavelength points will 
correlate. In other words, systematic errors will cause the 
neighbouring data points to tend to be on the same side 
of the calculated curve. Not surprisingly, the linear corre- 
lation coefficient r of neighbouring residuals, i.e. the sys- 
tematic error, was significantly larger for the solution with 
the (false) solid body rotation than for the one with the 
(correct) differential rotation (Table ||: 1st and 2nd line). 



The systematic error could also be seen in the slightly 
larger deviation d. 

To show that this criterion (lowest correlation of resid- 
uals) does not favour any particular spot distribution, we 
also simulated synthetic observations with a polar spot 
and no differential rotation (CASE II). Again the simu- 
lated data was used to retrieve temperature maps with 
and without differential rotation. The correct stellar pa- 
rameters gave a solution consistent with the input map. 
In the solution with artificial differential rotation the po- 
lar spot was reduced, as expected. But also in this case, 
an erroneous differential rotation coefficient caused corre- 
lation in the residuals and increased the deviation (Table 
|: 3rd and 4th line). 

Table 8. The quality of the surface imaging solution for 
rigid body (a = 0) and antisolar differential rotation (a = 
—0.2) in CASES I and II: mean deviation d of the solution 
from the synthetic spectroscopic "observations" and linear 
correlation r of subsequent residuals of the line profiles 



Input parameters 




Solution 






polar spot 


v sin i 


a 


v sin i a 


d 


r 


no (CASE I) 


70 


-0.2 


70 -0.2 


0.30 % 


0.035 


no (CASE I) 


70 


-0.2 


72 0.0 


0.33 % 


0.198 


yes (CASE II) 


72 


0.0 


70 -0.2 


0.32 % 


0.110 


yes (CASE II) 


72 


0.0 


72 0.0 


0.30 % 


-0.001 



We conclude, that aiming for a maximum correlation 
between images obtained from different spectral lines will 
favour solutions with polar spots. Furthermore, systemat- 
ical errors in the calculated lines profiles increase, if the 
value of a is erroneous. 

In the simulated tests we correlated the errors of neigh- 
bouring points. But in our real data of HD 199178 neigh- 
bouring residuals will anyway be correlated, because the 
spectra have been rebinned to a common wavelength scale. 
This effect can be avoided by comparing each point in the 
spectra to the point two steps ahead. The correlation of 
next-to-neighbouring residuals was checked for the July 
1995 surface imaging solutions with the differential coeffi- 
cients a = -0.17, a = 0.0 and a = 0.15 (Fig. |). 

The results in Table || showed that the solution with 
a = —0.17 gave both the best fit (lowest d) and smallest 
correlation of the residuals (lowest r). However, we note 
that the correlation was highly significant even for the 
best stellar parameters (usin?=70 km s _1 , a = —0.17, 
i = 60°). Apparently there still remained systematic errors 
in the solution. This could also be seen from the fact that 
the deviation was larger than the observational noise. The 
most likely cause for these systematic errors were inaccu- 
racies in the calculated line profiles (missing lines, errors 
in spectral parameters or element abundances, non-LTE 
effects, etc.). Even though all systematic errors could not 
be eliminated, we consider negative surface differential ro- 
tation more likely than the polar spot, because the best 
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quality of the surface imaging solution was achieved with 
a = —0.17. Still, we emphasize that the a-value was used 
to adjust the calculated profiles to fit the observations. 
Therefore our estimate for a may be biased by errors in 
the line calculations. 

The presence of differential rotation is supported by 
the variation of the photometric period, although pho- 
tometry alone can only give a lower limit for \a\. Inserting 
the means of the spot latitudes and photometric periods 
of Table |4| into Eq. ||, the value of a = —0.17 would give 
the rotation periods at the equator P G quator ~ 3 d 81 and 
the pole Ppoio ~ 3 d 25. The reliable pho tomet ric periods of 
HD 199178 determined by Jetsu et al. ( |l999b| ) were within 
the range 3 d 191 ± 0.017 < F phot < 3 d 



559 ± 0.017. Since 



HD 199178 is a rapid rotator, one would expecte d it t o 
have only high latitude spots (Schiissler & Solanki 1992| ). 
Therefore it is natural, that we do not find the equatorial 
rotation period from photometry. But the fact that our 
-Ppoic is larger than the smallest measured P p hot indicates 
a slight discrepancy. 

By combining the time series analysis of photometry 
with individual surface images we can study if the pho- 
tometric rotation period depends on the latitude of the 
main spot. Antisolar differential rotation would be seen 
as a shorter period for higher latitude spots. We note that 
this is indeed the case for the photometric periods and 
latitudes of the spot in Table [|. But since we have only 
three measurements, this tendency is far from statistically 
significant. Furthermore, the differences in the spot lati- 
tudes obtained from the surface images are too small to 
draw reliable conclusions. 

There are other methods to measure the surface dif- 
ferential rotation. The approach most commonly used in 
Doppler imaging, is to cross-correlate constant latitude 
slices of images and calculate the shift in longitude as a 
funct ion of the latitude (cf. Donati & Collier Cameron 
1997 ). This approach requires that there really are long- 
lived spots at different latitudes. The method may also 
fail if there are too many artifacts in the images. Artifacts 
will often appear at nearly the same longitude as the 
real features or as vertical stripes in the images. A cross- 
correlation analysis of such images will give a nearly rigid 
body rotation. In our images, there is only one major spot, 
and the low latitude features are probably artifacts. For 
this reason, we do not attempt to cross-correlate our im- 
ages. 

Antisolar differential rotation has previously been 
found in some active late- type stars. Johns-Krull ( 1996 ) 
analyzed the rotationally broadened spectral line profiles 
of three T Tauri stars and got the best fits with antisolar 
differential rotation. But he also noted that a polar spot 
could produce this effect. Weak antisolar differential ro- 
tation has been found in Doppler images of the RS CVn 



stars UX Ari (Vogt & Hatzes 1996), HR 1099 (Vogtet al 



199S) and HD 106225 (Strassmeier 1994: Hatzes 1998). 



7. Conclusions 

The surface images of HD 199178 presented here resem- 
ble the results rece ntly o bta ined f or FK Comae Berenices 
(Korhonen et al. |l999j & fcOOOQ . The dominating fea- 
ture is a long-lived high latitude spot (or spot group), 
with slightly changing position, shape and strength. Such 
changes could already be seen in the July 1994 and August 
1994 images. In the November 1994 image, the spot had 
weakened and partly disintegrated. In July 1995, the spot 
had recovered its strength. The changes in the photomet- 
ric light curves of HD 199178 during 1994 and 1995 can 
thus be explained by evolution of a single large spot rather 
than by abrupt changes of the whole spot configuration. 

The long-lived spot structure suggests there is a sta- 
ble active longitude in HD 199178. The high latitude of 
the spot could be explained by that the Coriolis force 
dominates over the buoyancy force in the magnetic flux 
tubes erupting from the deep parts of the convection zone 
(Schiissler & Solanki 199S ) . The changes in the spot lati- 
tude could then be connected to variations in the magnetic 
field strength. The changes in the latitude will in turn lead 
to drifts in the longitude if there is surface differential ro- 
tation. 

In previous temperature maps of HD 199178, the "flat 
bottomed" absorption line profiles have been interpreted 
as a signature of a large and very cool polar spot (Vogt 
1988; Strassmeier et al. 1999). We achieved the best fit 
to the spectral observations by introducing an antisolar 
surface differential rotation of a = —0.17. With a > 0, a 
polar spot appeared in the images. We considered the so- 
lution with antisolar differential rotation to be more likely 
than the polar spot, because it reduced both the discrep- 
ancy between the observed and calculated line profiles, as 
well as the correlation between residuals of neighbouring 
wavelength points (i.e. systematical errors in the calcu- 
lated line profiles). 

During 1994 and 1995 the photom etric rotation pe- 
riod varied by 3 %. Jetsu et al. ( 1999b ) interpreted these 
variations as a signature of surface differential rotation. 
Our images show that the photometric minimum of the 
light curves was apparently caused by the same persistent 
high latitude spot (group). Therefore the changes in the 
photometric period could not be explained by a changing 
pattern of spots. We also note that in our surface images 
the changes in the spot latitude were less than 10 degrees. 
The differential rotation would have to be strong in order 
to cause the observed changes in the photometric rotation 
period. 

There are certainly other alternative explanations for 
the flat bottomed line profile shape than antisolar differ- 
ential rotation. The very rapid rotation of FK Comae- 
type stars, perhaps being recently coalesced binaries, sug- 
gests there could be significa nt dev iations from a spherical 
shape (cf. Weltey & Ramsey 1994 ). Our understanding of 
absorption line formation in the photosphere of these ex- 
ceptional stars may also be insufficient. In particular, the 
strong absorption lines may be subject to non-LTE effects. 
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On the other hand, Bruls et al. (1998) showed that chro- 
mospheric activity and non-LTE effects could only explain 
the box- like shape of some (e.g. Fe I A 6430 A), but not 
all (e.g. Fe I A 6411 A and Ca I A 6439 A) of the lines 
used in our analysis. 

Weak antisolar differential rotation has previously 
been detected in some RS CVn stars (cf. Vogt et al. |l999| ). 
These detections were made by comparing Doppler images 
having a time difference of up to several months. This 
approach is, however, only possible for stars with a very 
small differential rotation, since the largest sensible phase 
difference is 0.5. Secondly, it is hard to confirm that the 
same particular spots are really present in different images 
taken months apart. 



Strassmeier et al. (1999) cross-correlated constant lati- 
tude slices of images to determine the differential rotation 
of HD 199178. They could not see any effects of differen- 
tial rotation but stated that this "non-detection does not 
necessarily mean there is no differential surface rotation" . 
We stress that the cross-correlation will provide reliable 
results only if there are long-lived spots on different lati- 
tudes. That is not the case in our images of HD 199178. 
Furthermore, the result may be biased by artifacts in the 
images. 

Differential rotation can also be detected by combining 
Doppler images and time series analysis of photometry. If 
there is a clear correlation between the latitude of the main 
spot and the photometric period, one could determine the 
law of the differential rotation. But this kind of analysis 
would require more data than was available for the present 
study. 
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